Double occupancy as a universal probe for antiferromagnetic correlations 
and entropy in cold fermions on optical lattices 
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We verify signatures of antiferromagnetic (AF) correlations in the double occupancy D [Gorelik et 
al., PRL 105, 065301 (2010)] and study their dimensional dependence using direct quantum Monte 
Carlo in dimensions d = 2,3 and Bethe Ansatz in d = 1. We find quantitative agreement with 
dynamical mean-field theory (DMFT) in the cubic case and qualitative agreement down to d = 1. 
As a function of entropy s — S/{NkB), D is nearly universal with respect to d; the minimum in D{s) 
approaches s « log(2) at strong coupling, as predicted by DMFT. Long-range order appears hardly 
relevant for the current search of AF signatures in cold fermions. Thus, experimentalists need not 
achieve s < log(2)/2 and should consider lower dimensions, for which the AF effects are larger. 
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A thorough understanding of materials with strong 
electronie correlations is not only desirable on intellee- 
tual grounds, but also due to their increasing technolog- 
ical importance, e.g^ in magnetoresistive and supercon- 
ducting devices |2[ . Theoretical investigations of cor- 
responding Hubbard type models using direct analytical 
techniques, numerical approaches for finite clusters, and 
the dynamical mean-field theory (DMFT) have already 
shed light on many strong-coupling phenomena includ- 
ing metal-insulator transitions, heavy-fermion and non- 
Fermi-liquid behavior, and various types of magnetic and 
orbital order Q. However, there are still important open 
questions, most notably regarding high-temperature su- 
perconductivity for which so far no mechanism could con- 
clusively be established. In this situation, the recent ad- 
vent of a novel class of correlated Fermi systems, namely 
ultracold fermionic atoms (such as ""^K and ^Li) on op- 
tical lattices, have opened a new promising direction of 
research: cold atoms are predicted to serve as quantum 
simulators for the Hubbard type solid-state Hamiltonians 
of interest 0, H]- 

Indeed, within a few years after the first achievement 
of quantum degeneracy in (single flavor) fermionic atoms 
on optical lattices , the Mott metal-insulator transition 
(MIT) was observed in two-flavor mixtures, based on sig- 
natures in the compressibility Ml and a suppression of the 
integrated double occupancy j8|- As a result, it is now 
established that the single-band Hubbard model 
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(with hopping amplitude t, onsite interaction U, and 
fiia = c\^c^^) can be realized to a reasonable accuracy 
using ultracold fermions in the interesting interaction 
range, which certainly supports the hopes of accessing 
also less understood Hubbard physics in similar ways. 



However, all attempts of realizing and detecting quan- 
tum magnetism in cold lattice fermions have failed so far. 
In fact, it has not even been possible yet to verify specific 
signatures of antiferromagnetic (AF) correlations which 
are ubiquitous in correlated electrons and believed to play 
an important role in high-temperature superconductivity. 
This type of physics clearly has to be under control be- 
fore cold fermions can really play a useful role as quantum 
simulators. Up to now the failures to detect AF signals 
have primarily been attributed to cooling issues [l^l . 
Indeed, the coldest systems achieved so far have central 
entropies per particle of s = S/{NkB) ~ log(2) while 
AF long-range order (LRO) on a cubic lattice is expected 
only for s < log(2)/2 [l3,[l3,[ll. 

In this Letter, we will argue that this discrepancy is not 
really relevant for the experiments currently performed or 
prepared in this context: both modulation spectroscopy 
T^ l and the superlattice approach 
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nearest- neighbor (NN) spin correlation function {&i ■ arj) 
(for Pauli matrices a). This is also true for the double 
occupancy D = {hi^hii) at large U/t (and temperature 
T = 0)M- 
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However, as we will show using the example of D, all such 
obscrvables are too local to be sensitive to LRO; given 
typical signal to noise ratios it seems very unlikely that 
the Neel transition could be detected in cold atoms in 
this way (assuming low enough s is finally reached). As 
a consequence, full dimensionality (i.e. an isotropic cubic 
optical lattice) is not essential; in fact, a restriction of the 
atoms to planes or chains might enhance AF signals. 

In the following, we will first briefiy recall the DMFT 
scenario put forward in and discuss arguments 

Hil l against the reliability of DMFT in dimension 
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FIG. 1. (Color online) Top: Double occupancy D{T) as es- 
timated from DMFT (circles), QMC (diamonds), and DCA 
[l^ (dashed lines) for (a) cubic and (b) square lattice. Thin 
hues: nonmagnetic DMFT results for T < T^^^'^ (arrows). 
Bottom: (c) The Neel temperature Tn{U) (determined by 
QMC [131 and DCA ^) is up to 30% below the DMFT 
estimate for d = 3. (d) While Tn vanishes in d = 2, the 
spin-crossover temperature [2^ agrees with T°^^^(i7). 



d = 3. Wc will show, by comparisons with direct dc- 
terniinantal quantum Monte Carlo (QMC) simulations 
[iol l , that the AF signatures predicted by DMFT survive 
even on the square lattice (d = 2) and arc surprisingly 
precise, up to rounding effects, in the cubic case (d = 3). 
Finally, we will focus on a coupling strength correspond- 
ing to the ground state Mott transition l2l| ([/ /t = 15 
for the cubic lattice) and demonstrate the effects of vary- 
ing dimensionality from DMFT (exact for d oo [23[) 
via d = 3 and d = 2 (from QMC) to d = 1 (from a 
thermodynamic Bethe ansatz, BA j24j). 

AF signatures in the double occupancy - A recent real- 
space DMFT study showed that the low-temperature 
formation of an AF core in a fermionic cloud on an optical 
lattice (with half filling = 1 in the center) is signaled, 
at strong coupling, by a significant enhancement of D in 
the same region. This DMFT scenario [l3| is reproduced 
for ri = 1 in Fig. [T^: at strong coupling U/t = 18 the 
double occupancy D{T) is asymptotically flat at low T 
in the nonmagnetic phase, but is strongly enhanced, by 
up to 75%, when AF order sets in below T^^'^'^ « 0.35t 
(arrow). The relative enhancement quickly decreases at 
smaller U and is lost at U/t « 10. The absolute en- 
hancement of D is largest for U/t 12 and should be de- 
tectable experimentally even in measurements integrat- 
ing over the whole cloud [l7| . 



This scenario was challenged recently [l8[ on the ba- 
sis of the dynamical cluster approximation (DCA) which 
relaxes the DMFT assumption of a momentum indepen- 
dent self-energy [2^: the DCA estimates of D (dashed 
lines in Fig. [T^) showed no clear AF related enhance- 
ment [l3l; however, these calculations could not enter 
the low-T AF phase. The reliability of DMFT estimates 
for D and s (in the nonmagnetic phase) at low T was also 
questioned based on comparisons with high-temperature 
expansions (HTE) [H] . It is indeed clear that the DMFT 
scenario cannot be correct in all aspects: after all, it is 
well-known that DMFT overestimates the Neel temper- 
ature by up to 30% in 3 dimensions (see Fig. [TJ;). Thus, 
the kinks in D{T) at T^^^^'^ are certainly unphysical - 
but is the whole scenario just a DMFT artifact? 

Comparison in d = 2 - For a first answer, let us turn to 
the square lattice (d = 2) for which the DMFT is a priori 
much less reliable than in d = 3. In fact, DMFT predicts 
AF LRO even in this case, with a maximum in T^^^"^ 
of about OAt at U/t « 8 (circles in Fig. [TJi), while the 
Mermin- Wagner theorem excludes LRO for T > Tjy = 0. 
However, in this case it is relatively easy to check the 
DMFT predictions (circles in Fig. [TJa) by direct QMC 
simulations of finite clusters, here of size 10 x 10 
(diamonds). After employing an (approximate) correc- 
tion for Trotter errors and verifying that finite-size ef- 
fects are negligible we consider this data essentially ex- 
act. The previously established accuracy of the DMFT 
at high temperatures H, 1^ evidently survives in d = 2, 
with no significant deviations from QMC for T/t > 0.8. 
Surprisingly good agreement is also found at low tem- 
peratures T/t < 0.2, although the stable DMFT solu- 
tions (circles) here correspond to the AF phase which 
at first sight appears unphysical. In contrast, DMFT 
calculations constrained to the nonmagnetic phase (thin 
lines) predict low-T features of D which are far off from 
the exact QMC data. This teaches an important les- 
son, relevant also for d = 3: paramagnetic phases which 
include short range order can be much more similar to 
AF phases (with AF LRO) than to nonmagnetic solu- 
tions (without any AF correlations). Of course, DMFT 
is still not perfect: QMC shows significant corrections of 
DMFT predictions, namely a rounding of the unphysical 
kinks, at T w T^^^'^. However, even T^^^^ has phys- 



ical significance: it (nearly) matches the spin coherence 
temperature (diamonds in Fig. [TJi) (ist - 

Comparison in d = 3 - QMC results on a cubic 
lattice, the target system of current AF related experi- 
ments, were obtained for clusters with 6"^ and 8'^ sites and 
carefully extrapolated to vanishing Trotter discretization 
At (large diamonds in Fig. [2]). These data show 
perfect agreement with the DMFT estimates (circles) at 
U/t = 12 both for T/t > 0.7 and for T/t < 0.4; thus 
the DMFT prediction of the D enhancement ~jl7l| is even 
quantitatively correct. Only at T^^^^ « 0.45t the QMC 
results smooth out the DMFT kink. 
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FIG. 2. (Color online) Cubic lattice (d = 3): Comparison 
of DMFT (circles) and direct QMC (diamonds) results for 
the temperature dependence of the double occupancy D{T) 
for the cubic lattice. DCA results taken from [T8|. Small 
diamonds: QMC results with finite Trotter discretization bias. 



At relatively weak coupling U/t < 8, signatures appear 
in the DMFT data in Fig. [2] which differ fundamentally 
from the strong-coupling scenario discussed so far: D{T) 
shows a broad minimum at T ss the rise towards lower 
T breaks down quite abruptly below T^^^"^ , remarkably 
approaching exponential fits to the high-T behavior (dot- 
ted lines). Apparently the system behaves as a (bad) in- 
sulator for T >t; the Fermi liquid behavior setting in for 
T < t enhances [13, but is destroyed below T^^'^^ by 
AF correlations. Also this weak-coupling DMFT scenario 
for D{T) is confirmed: QMC predicts (large diamonds) 



a peak right at T^'^^'^ « 0.2t for U/t = 4 and quickly 
converges towards DMFT for lower T. The deviations in 



the range Tj 



DMFT 



< T < t can be traced to developing 



AF correlations which already reduce the Fermi liquid 
enhancement of D. Note that (at U/t = 4) the discrep- 
ancies between QMC and DMFT arc much smaller than 
typical QMC discretization errors (data for Art =1/8: 
small diamonds) and that DCA (dashed line) apparently 
misses the AF physics at T/t < 0.2. 

LRO leaves traces in the QMC estimates of D{T) only 
in the weak-coupling regime U/t <6 where Tn « T^^^"^ 
(cf. Fig. [It). At strong coupling U/t = 12, D{T) does not 
show visible features at Tn « OM, which suggests that 
local spin correlations (which determine D and current 
AF observables [l^jllEl) are hardly sensitive to LRO and, 
consequently, dimensionality in this regime. 

Impact of dimensionality and entropy ~ In order to 
gain more insight into these issues, DMFT and QMC 
data for the cubic lattice are compared at [//t = 15 with 
QMC results for the square lattice and BA solutions of 
the infinite chain at (nearly) equivalent (29| interactions 
in Fig.[3K- Here, the DMFT data (circles) can also be in- 
terpreted as an exact result in infinite dimensions. After 
rescaling (29j . we find rapid convergence with increasing 
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FIG. 3. (Color online) Hypercubic lattice (1 < d < 3) at 
strong coupling: a) D(T) as estimated from DMFT (d = 3, 
circles), QMC (d = 2 3, diamonds), and BA (d = 1, 

dash-dotted line), b) Corresponding estimates of entropy per 
particle s — S/N. All interactions correspond approximately 
to the ground state Mott transition at U/{\^t) ~ 6. 



dimensionality at high T and generally similar shapes [30| 
of D{T) for 1 < d < oo. However, c? ^ 3 would appar- 
ently be needed in order to converge to the DMFT results 
also at T]^^^"^. Furthermore, the minimum in D{T) oc- 
curs at about twice T^^^^ in dimensions 1 < d < 3, 
reinforcing doubts about the usefulness of DMFT esti- 
mates of D for thermometry [l8| . 

It is well-known that nonmagnetic DMFT yields an 

entropy s '^~^°> log(2) (dotted line in Fig. [5)d), which is 
clearly unphysical [19[. However, the AF DMFT solution 
(circles for T < T^^'^'^) recovers the QMC results for the 
cubic lattice (diamonds) at T < Tn ~ 0.25t; remarkably 
the latter coincide with the Heisenberg limit of s(T) for 
T < Q.8t. In general, the dimensional dependence of s{T) 
nearly mirrors that of D{T). Thus, dimensional effects 
and DMFT errors should be minimal when using s as 
a (dimensionless) measure of temperature (which is of 
primary interest to experimentalists anyway). 

Indeed, as seen in Fig. IH D{s) looks strikingly similar 
in all dimensions; in particular, the minimum in D (at 
strong coupling) corresponds to s « log(2) = s^'^^'^ in 
all cases! While it might appear surprising that this be- 
havior persists down to d = 1 it is clear that s < log(2) is 
only possible for a two-flavor system at n = 1 by spin co- 
herence, i.e. the development of (possibly short ranged) 
AF correlations; these, in turn, enhance D (l7| . 

Thus, the evolution of D is a near-perfect thermome- 
ter for ultracold atoms measuring AF correlations and 
(in d = 3) the proximity to AF LRO. In fact, we would 
argue that any positive deviation of D{s) from the non- 
magnetic background (shaded in Fig. [4|) originates from 
AF correlations which are strong for s < log(2) w 0.7 in 
all dimensions (and coincide with LRO in c? = cxd and, at 
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FIG. 4. (Color online) Hypercubic lattice (1 < d < 3) at 
strong coupling: Double occupancy as a function of entropy 
per particle. In all cases, the minimum of the double occu- 
pancy corresponds to s ~ log(2) (dotted line). The shaded 
area indicates the nonmagnetic contribution to D. 

s < stv ~ log(2)/2, in d = 3) and are strictly zero only 
in the limit d oo for s > log(2). This enhancement is 
larger in lower d, consistent with the known dimensional 
dependence of (cri-a-j) for the Hcisenberg model at T ~ 0: 
(a, ■ &j)o = -1.00 Cd = oo, Weiss MF)^1.20 (d = 3) 
[3l|; -1.34 {d = 2) [33]; -1.77 (d = 1) [3]. Thus, irre- 
spective of the measurement technique, signatures of AF 
correlations may be easier to detect experimentally (at 
fixed s) for lower (effective) dimensionality. Conversely, 
a tuning of the hopping amplitude in z direction could 
help to discriminate magnetic effects from those of charge 
excitations; similar ideas including frustration will be ex- 
plored in a separate publication |33| . 

Conclusion - In this Letter, we have demonstrated 
that DMFT predicts temperature dependencies of local 
quantities, e.g. D{T), more accurately in d = 3, d = 2 
than expected. Especially, both the AF induced enhance- 
ment (at strong coupling), and the suppression (at weak 
coupling) of D survive. The temperature scale given by 
jiDMFT^jj-j corresponds to a spin-crossover in finite di- 
mensions. As a function of entropy per particle s, the 
double occupancy is nearly universal with respect to di- 
mensionality; in particular, the minimum in D{s) always 
occurs at s w log(2) at strong coupling, as predicted by 
DMFT. Thus, we have established a prominent and spe- 
cific signal of AF correlations in an entropy range that is 
in immediate experimental reach, with the prospect of ex- 
tending the use of D(s) for thermometry [11| to the most 
interesting range s < log (2). Our results also validate 
the RDMFT approach [3J| for quantitative simulations 
of inhomogeneous 3-dimensional systems. 
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